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We analyse the dynamics of spinodal decomposition in inflationary cosmology using the closed 
time path formalism of out of equilibrium quantum field theory combined with the non-perturbative 
Hartree approximation. In addition to a general analysis, we compute the detailed evolution of two 
inflationary models of particular importance: A3> 4 new inflation and natural inflation. We compute 
the metric fluctuations resulting from inflationary phase transitions in the slow roll approximation, 
showing that there exists a regime for which quantum fluctuations of the inflaton field result in 
a significant deviation in the predictions of the spectrum of primordial density perturbations from 
standard results. We provide case examples for which a blue tilt to the power spectrum (i.e. n s > 1) 
results from the evolution of a single inflaton field, and demonstrate that field fluctuations may result 
in a scalar amplitude of fluctuations significantly below standard predictions, resulting in a slight 
alleviation of the inflationary fine tuning problem. We show explicitly that the metric perturbation 
^ ■ spectrum resulting from inflation depends upon the state at the outset of the inflationary phase. 
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(S| 1 I. INTRODUCTION 
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In recent years, there has been a strong interest in the dynamics of quantum fields in the early universe. This 
interest has led to a better understanding of a number of processes including the formation of topological defects 
during early phase transitions jjj , the reheating of the universe after inflation Q , and the dynamics of inflation itself 
H . In the particular case of inflationary reheating, our improved understanding has been revolutionary and has 
"sj" ' significantly reshaped the subject Q. 

The lessons provided by these studies are varied. One crucial aspect is the importance of using time-dependent 
techniques to study processes of the early universe. It has been repeatedly shown that classical and one-loop effective 
potentials are poorly defined and of little use in dynamical systems; they should only be used to determine static 
quantities such as the ground state of the system ||. Another common theme is the importance of non-linear 
' corrections to the linear dynamics. These corrections have been found to be quite dramatic in studies of phase 
transitions and reheating. 

Despite these important advances in studies of quantum fields in the early universe, it is still widely believed that 
these techniques have little to add to our understanding of the inflationary phase itself. The belief that the inflaton 
• follows a classical trajectory determined from the classical effective potential with only perturbatively small quantum 
corrections [|| is still widely held. While the techniques of stochastic inflation allow these corrections to add up, 
they do so in an incoherent fashion through a repeated summation of one-loop effects without self-consistent inclusion 
of higher order corrections Q . 

Much work has been done to verify that the dynamics of the inflaton field is predominantly classical. The existence 
of a particle horizon and the natural squeezing of states due to the near exponential redshifting of field modes justifies 
a quasi-classical description of the inflaton field up to a perturbatively small component of field fluctuations with 
wavelengths inside the de Sitter horizon However, it is important to emphasize that the validity of a classical 
evolution of the inflaton field alone does not justify the use of the classical effective potential for a non-linear dynamical 
system. As has been shown quite clearly by the classical field theory simulations of the early stages of reheating, the 
full non-linear dynamics of even a purely classical system depends strongly on coherent effects of backreaction due to 
the field's fluctuations JTo[ ] . It is then not unreasonable to expect that non-linear field fluctuations during inflation 
could result in a departure from the dynamics derived directly from the classical effective potential. 

In this report we address these issues using the techniques of out-of-equilibrium quantum field theory. We examine 
a class of models in which strongly quenched inflaton evolves under the influence of a negative mass squared in the 
potential, a process which, following the terminology of such phase transitions in condensed matter physics, we refer 
to as spinodal decomposition p]]. This class of inflation models includes new [jl2| and natural inflation |Q, as well 
as many models of hybrid inflation jlJJ . Such inflation models are of particular interest in the present study because 
the evolution from the initial to the final state of the system is necessarily a non-linear process. 

The particular case in which the inflaton field is treated as a component of an O(N) vector in the large N limit 
has already been detailed for the case of new inflation, where it was found that the full non-perturbative quantum 
dynamics does in fact reproduce an effectively classical trajectory for the evolution of the inflaton. The growth of 
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quantum fluctuations results in a dynamical flattening of the potential |15| , an analogue of the Maxwell construction 
commonly used in studies of the equilibrium properties of phase transitions. In this earlier work, it was found that 
the effectively classical trajectory of the inflaton field, a result of the phenomenon known as zero mode reconstruction, 
is a consequence of Goldstone's theorem fl6f| . 

However, the case of a single, real inflaton field is quite different. Here, the long wavelength quantum fluctuations 
of the field also reassemble themselves into a semi-classical field, but in this case the assembled field does not obey the 
classical equations of motion of the original potential. Rather, the inflaton may be broken up into two components. 
The first is the mean field <p and obeys the classical equation of motion expected from the original potential except 
that it is coupled to a second field. This second field cr, constructed through the assembly of quantum fluctuations, 
obeys a modified equation of motion. The result is an effectively classical theory of two coupled fields, refered to as 
spinodal inflation |l7| . 

The observational consequences are dramatic. As there are effectively two fields, the evolution becomes quite 
complicated and depends on the initial conditions. There are two regimes. In the first, the mean inflaton field, defined 
as the expectation value of the quantum field, has an initial value greater than the expansion rate <j>(to) > Hq/2it and 
the semi-classical field a never becomes dynamically relevant. The evolution reproduces all of the standard results 
for that particular model of inflation and we can think of this as the classical regime. In the second, quantum, regime 
for which 0(£o) < Hq/2tt, the influence of the a field is quite important. In this case, observational quantities such as 
the amplitude and spectrum of primordial density fluctuations depend not only on the parameters of the model, but 
also on the particular value of the initial mean field 4>{to). 

A simple single field model can therefore produce a range of observational results for any given choice of parameters. 
In fact, due to the effective two field dynamics, it is possible to produce observational features not possible in the 
simple classical, single field version of the same theory, such as the generation of a blue primordial power spectrum; 
this can occur much in the same way as it does in hybrid inflation. 

We begin with a short introduction to spinodal models of inflation and the need for a fully out-of-equilibrium and 
non-perturbative description of the dynamics. We write down the self-consistent Hartree equations of motion for a 
general spinodal potential, followed by an explanation of the assembly of quantum fluctuations and how this results 
in an effectively classical two field model. Next, we move on to a detailed analysis of the two most important spinodal 
models of inflation, A^ 4 new inflation and natural inflation. 

Having determined the evolution of the field, we wish to examine the observational consequences of spinodal 
inflation. Beginning with the gauge invariant formulation of gravitational perturbations of Mukhanov, Feldman, and 
Brandenberger we note that while the effective dynamics is that of two fields, returning to the full quantum 
theory we see that only one set of field fluctuations couples to gravity perturbations. 

We then provide a complete and detailed numerical analysis of the dynamics of spinodal inflation in single field 
models, including computations of the primordial spectrum of scalar and tensor perturbations which result. The 
consequent C; spectra for a particularly interesting example is shown to allow for the direct relation of spinodal 
effects to observation. Here we see explicitly the exciting features possible through the dynamics of spinodal inflation. 



II. SPINODAL MODELS OF INFLATION 

We envision the evolution of a scalar field near the top of a potential of the form sketched in Fig. |l|. The potential 
may be expanded in the form 

V(*)=K-\n 2 & + ±& + -. , (2.1) 

where the cosmological constant contribution K is chosen such that the potential is zero in the true vacuum, and /i 2 
is positive. 

Initially, the field will "roll" slowly toward one of the minima of the potential. This is the regime in which inflation 
will take place. To a first approximation, we ignore the quartic term and see that the initial evolution follows that 
of a free field in an inverted harmonic potential. This evolution has been studied in great detail in the context of 
inflationary cosmology |20). For early times, the field grows exponentially. Eventually the higher order terms in the 
potential become important, with the result that any perturbative analysis of the dynamics will break down and must 
be augmented by some non-perturbative technique. 

Our choice of approximation is further restricted by the fact that the system we wish to study is not in thermal 
equilibrium, thus leading us to real-time methods. We emphasize that equilibrium constructs, such as the effective 
potential, are completely inadequate tools for this problem. 
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FIG. 1. A typical potential with a region of negative curvature. The dashed line is the spinodal line, separating the unphysical 
spinodal region (above the line) from the physical re gion (at and below the line). The exact shape of the dotted portion of the 



curves depends upon the higher order terms in Eq. (2.1) 



The simplest approximation satisfying the requirements is the Gaussian variational approximation, in which the 
quantum density matrix is restricted to take on a Gaussian form. Also known as the Time Dependent Hartree-Fock 
approximation, such mean field techniques have been utilized in quantum mechanics dating back to Dirac jgjj . It is a 
standard technique in chemistry, condensed matter physics, and nuclear physics and has led to a better understanding 
of the structure of a number of phase transitions. 

A. Real Scalar Field: Hartree Dynamics 

In what follows, we assume a spatially flat Robertson- Walker metric: 

ds 2 = dt 2 - a 2 (t)dx 2 . (2.2) 
We now derive the equations of motion for a real scalar field with Lagrangian 

L = lv^(x)V^(x) - V [*(x);t] , (2.3) 



within the self-consistent Hartree approximation |22|. We break up the field $ into its expectation value plus a 
fluctuation about this value: 

$(x,t) = <j>{t)+^{x,t) , (2.4) 

0(t) = mS,t)). (2.5) 



Here, <j) depends only on time due to space translation invariance as is consistent with the metric (2.2). By definition 
{iP(x,t))=0. 

The Hartree approximation consists of replacing ip 2n by ci(ip 2 ) n 1 tp 2 +C2(ip 2 ) n and ip 2n+1 by C3(ip 2 ) n ip, where the q 
are constant factors whose values are determined by Wick's Theorem. This Hartree factorization may be summarized 
as follows: 

^ _j2n)I { 2 x a _ (2n)!(n-l) a 

ij2n+1 ^ (2n+l)! 

2 n n\ 

Given this factorization, any function F(<&) becomes 
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where we use the notation 
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F{cj>). 



(2.7) 



(2.8) 



Note that the latter two terms on the right hand side of Eq. (2.7) have zero expectation value. We therefore find that 
the expectation value of a function factorizes as 



71=0 ' ^ ' 

The equations of motion for the mean field 4> are given by the tadpole condition (ijj) = 

*) + 3f + f;i(<|>)V<-»(^o, 



(2.9) 



(2.10) 



where we have used the metric (2.2) and have factorized the expectation value (V ($)) according to Eq. (2.9). We 
define the Fourier transform of the Wightman function by the expression 



G{x, t;x ,t) 



d 3 k 



e <&(*-*)G fc (t,t') , 



where we have used the property of space translation invariance. The quantity {ij) 2 (i)) = —iG(x,x) is constructed 
from the mode functions obeying the equation 



[G k (t,t)]- 1 f k {t) = Q 



(2-11) 



together with the appropriate Closed Time Path boundary conditions 0]. The operator [G k ] 1 is given by the 
quadratic form appearing in the generating functional. Explicitly, the f k (t) obey 



d 2 a(t) d k 2 
dl? + ^(tjdi + ^ 2 {t) 



oo 1 

T- 

71—0 
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/*(*) = o 



(2.12) 



where we have again used the factorization ( [2.9] ) to express the potential term (V ($))• 

As mentioned above, the quantity (ip 2 ) is determined from the mode functions f k combined with Closed Time Path 
boundary conditions appropriate to the chosen initial state. For an initial state in thermal equilibrium at an initial 
temperature T, we have 



d 3 k 



l/ fc W| 2 coth(^|) 



(2.13) 



Note that in the zero temperature vacuum state given by the T — » limit, the hyperbolic cotangent has the value 1. 
The frequency to k appearing here is given by 



k 2 + a 2 {ta) ylfm ))Y v j M) 

i n V 2 / 

n=0 v 7 



mo)) 



(2.14) 



For the case of the A3> 4 theory, there is an additional term proportional to the Ricci scalar — a 2 (to)R(to)/6 which 
appears on the right hand side of this expression for uj k . This term arises when one considers initial conditions 
corresponding to the adiabatic vacuum state in conformal time pq , pq| , which is necessary if we wish our initial 
vacuum state to match the Minkowski vacuum in the limit a(t) — > 1. However, this term is not necessary if the scalar 
field is taken to be part of a low energy effective theory, such as is the case with the natural inflation case we analyze 
below. 

Using the vacuum state for the mode functions defined by the initial frequency spectrum of Eq. (2.14) leads to the 
following initial conditions on the f k (t): 
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fk(t ) = 



/,(/„■)= (-^j-iwfcJ.A (V). 



(2.15) 



A final note is tha t the initial frequencies given by Eq. (2.14) may be imaginary for low k modes. In this case the 
initial conditions ( 2.1 5| ) need to be modified for low k. This may be done in a variety of ways with little effect on 
results [M. Here, we choose a smooth interpolation between low k modes with modified frequencies and the high k 
modes which remain in the conformal vacuum state with frequencies u>k ( |2.14 ) : 



M tanh 



■M' 



\M 2 



where 



n=0 ' ^ ' 



(2.16) 



(2.17) 



This completes our set of equations of motion of the matter fields within the Hartree approximation. 



B. Gravitational Dynamics 



We will treat gravity in the semi-classical approximation, in which the expectation value of the full quantum energy- 
momentum tensor acts as a classical source to the Einstein gravitational tensor. The semi-classical Einstein's equation 
reads 
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(higherderivatives) = — (T^ 1 ) , 



(2.18) 



where Gjv = 1/Mp t is (bare) Newton's constant, A is the (bare) cosmological constant and the components of the 
Einstein curvature tensor using the metric (2.1) are: 



-R 



M a - + a - 

a 



(2.19) 
(2.20) 



The higher derivative terms included in Eq.(2.18) are needed for renormalization purposes; Newton's constant and 
the cosmological constant will also be renormalizcd (see below). 

Once again, we use the factorization (2.9) to determine the right hand side of (2.18). Defining the additional 
integrals 



<(v^)) 2 ) = 



d 3 k 
(2tt)" 3 
d 3 k 



-|/ fe (t)| 2 coth 



(-) 

\2TJ 



k 2 \ f k (t)\ 2 coth 



f<£k \ 
\2T) ' 



we find for the energy density e and the trace T of the energy momentum tensor: 

t = <r;> = -tf - {^) + i ((v^) 2 ) + 4 f; i ( ) n v^\<t>) . 

a T l = U - \ / 



(2.21) 
(2.22) 

(2.23) 
(2.24) 



The pressure p is arrived at from these expressions through the relation p = (e — T)/3. The equation of state of the 
system is characterized by the quantity w — p/e. 
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C. Regularization and Renormalization 



The mode integrals appearing in Eq. ( 2.13 ), ( 2.21 ). and ( [2.22 ) are formally divergent and must be regularized in 
order to perform any practical computation. 

In the special case of a renormalizable potential V($>), for example in the A$ 4 model, we would like to fully 
renormalize the theory. Our choice of renormalization procedures is somewhat limited by the requirement that the 
dynamics be amenable to numerical analysis. However, a number of groups have recently addressed this problem 
either by means of adiabatic regularization with a simple cutoff at large momentum [27 25 2£|, as was first developed 
by Paul Anderson p{| , or by using a scheme based on dimensional regularization J3C |7~^ 

In practice, we use the simple scheme developed by the Pittsburgh-Paris collaboration. This scheme has the 
advantage of being very easy to implement and it has the attractive feature that all subtractions are absorbed into 
counterterms renormalizing the bare couplings in the equation of motion and the Friedmann equation (the latter 
must be extended to include a cosmological constant and a higher order curvature term). However, we mention that 
it does not include the finite subtractions which would be necessary to give the correct conformal anomaly. These 
terms, the finite subtractions detailed in p7| and [ p8| , are formally important, but in the present context, as the 
inflaton self-coupling is typically required to be of order 10 -12 or smaller, such terms will have absolutely no influence 
on numerical simulations as their contributions are much smaller than the numerical accuracy of the computations. 
In fact, although we do not do so here, it is normally safe to drop the logarithmically divergent terms from the 
simulations as well without influencing the results. The simulations are checked after the fact to verify that they 
satisfy the covariant conservation of the energy-momentum to within their numerical accuracy, and to ensure that the 
results are independent of the value of the momentum cutoff. 

For models with non-renormalizable potentials, such as natural inflation, we have to satisfy ourselves with the treat- 
ment of the model as a low energy effective theory with a well defined cutoff. Again, we implement the regularization 
by means of a large momentum cutoff. 



III. EARLY TIME DYNAMICS AND REASSEMBLY 



As mentioned above, the initial linear dynamics in spinodal models of inflation is well understood. This period is 
characterized by exponential growth of both the mean field (f> and those mode functions fk with physical wavelength 
greater than the Hubble distance, that is, k/a(t) H(t) = a(t)/a(t). Since H(t) is approximately constant during 
inflation while a(t) is growing exponentially, clearly more and more modes satisfy this condition at each subsequent 
time. However, the very long waveleng th m odes which crossed the horizon very early on in the evolution will tend 
to dominate the quantity (ip 2 ) of Eq. ( 2.13 ) simply because they have experienced the spinodal instability for the 
longest time. This is a very important characteristic of spinodal inflation which sets it apart from other models: the 
dynamics is driven by a super-horizon scale quasi-particle condensate. On scales smaller than the horizon, it is not 
possible to distinguish such a condensate from a purely homogeneous background field. Any possible measurement 
will determine only the combined properties of the condensate and the mean field. 

Given the assumption of initial conditions near the local maximum of the potential, and provided with a finite 
renormalized or regularized two-point function along with very small values for the higher order couplings in the 
Lagrangian, the early dynamics is well approximated by a linear analysis. The equation for the mean field for any 
potential of the form (12.11) is simply 



4>{t) + ZH{t)(j>{i) - ii 2 cj){t) = 



(3.1) 



To this order, we may take the Hubble parameter to be constant, H(t) = Ho, in which case the solutions to this 
equation are simple exponentials. Only the growing term is relevant, so that we have the early time solution 



4>{t) 



Do exp 




(3.2) 



The value of cf>o depends on the precise initial conditions for ip and 
The mode functions obey the similar equation 



f k (t) + 3H(t)f k (t) 



(3.3) 



which for constant H(t) = Hq and corresponding exponential a(t) = exp(i?o£) has the solutions 
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f k {t) « exp 



;H t 



-H t 



BkJ-v I -r—e 
Ho 



-Hot 



(3.4) 



The constants Ak and Bk are determined by the initial conditions on the mode functions. These solutions oscillate 
with an envelope proportional to l/a for sub- horizon modes (with k/a > Ho), but the solutions are growing and 
decaying exponentials in the opposite limit k/a <C Ho, leading to the statement that super- horizon modes have 
an exponential instability. In this limit, we may again discard the exponentially decaying term. We then have for 
super-horizon modes 



C, 



exp 



(3.5) 



Here, Ck is roughly equal to the value of the mode fk evaluated at the time t(k) that the mode crosses the horizon as 
determined by the condition k = cxp[Hot(k)]Ho- For modes initially far within the horizon the general dependence of 
Ck on k may be estimated from the initial behavior of the mode functions fk ~ k~ x / 2 and the decay of the envelope 
of th e Bessel solutions which provides the standard result that Ck oc fc~ 3 / 2 . We note tha t wh ile the exponential form 
(3.5) is only an asymptotic solution for small arguments of the Bessel functions of Eq. (3.4), due to the exponential 



behavior of this argument it very accurately describes the evolution any mode function within a Hubble time of 
horizon crossing. 

We are now in position to compute an expression for the condensate (tp 2 ). By separating the momentum space 
integral over super- and sub-horizon modes respectively, we can take advantage of the expression ( pl| ) to find for 
early times 



(?A 2 ) = exp 



(,H 



d 3 k \c k \ 2 , ( ^ 2; 



2(2tt) 3 k 



3 u2u-3 



k>aH 



(3.6) 



The latter, sub-horizon term contains all subtractions due to renormalization. After a few Hubble times, it is safe to 
neglect this term compared to the exponentially growing super-horizon contribution. To determine which modes are 
most important to the evolution, it is convenient to examine the contribution of each squared mode per logarithmic 
momentum interval, dk/k. Using the k~ 3 ^ 2 behavior of Ck, we find for modes which are originally far inside the 
horizon, but which have since crossed outside, that this contribution is proportional to k~( 2,J ~ 3 \ Since 2v > 3, the 
integral is dominated by those modes which crossed the horizon the earliest. This is important for any numerical 
analysis as it allows one in practice to set a cutoff in the calculation of (ifj 2 ) with well controlled errors, avoiding the 
problem of including a number of mode functions which grows exponentially with the scale factor. 

As a consequence of the formation of the condensate, it becomes possible to form an accurate and simple model of 
the complete system, in which the full two-point fluctuation (-0 2 ) is replaced by a nearly homogeneous, and effectively 
classical field. This produces a model in which two effectively homogeneous classical fields, the mean field coupled to 
a fluctuation field, accurately describe the dynamics. The condensate field is defined as 



a(t) = V(^(t))k<aH 



(3.7) 



For early times, it is given by the first term on the right hand side of Eq. (3.6). The square root of the value of the first 
integral of (3.6) a few Hubble times after the beginning of inflation provides an effective initial condition on a. For a 
zero temperature initial state, which for simplicity we take to be the case in what follows, this is found numerically 
to be of order H/2tt pTq ] . The effect of a finite temperature initial state is to increase this value by a factor of y/T/fi 
without modifying any of the qualitative features described in this study. The dynamics modelled by this classical 
field will be accurate up to perturbatively small corrections to the dynamics due to the sub-horizon modes contained 
in the final term of Eq. ( |3.6| ) . 

This system of effective homogeneous fields is referred to as a reassembled system. We will refer to the <fi field, 
the expectation value of the full quantum field, as the mean field, while we refer to the semi-classical a field as the 
fluctuation or condensate field. 

While this condensate forms during the linear regime, eventually the dynamics becomes non-linear. It is this non- 
linear evolution which we wish to study. We will be particularly interested in examining how the interaction of the 
condensate with the mean field can influence the dynamics. 

We consider primarily the dynamics of single field models (the large N case was studied previously |l(Jj). Here, 
the interactions of the condensate and the mean field will be seen to lead to a complicated evolution in which initial 
conditions play a primary role. 
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IV. NON-LINEAR DYNAMICS 



We begin with an analysis of the stationary solutions for the mean field, which obeys the equation of motion (2.10). 
There are two primary late time solutions. The first is the trivial solution with <fi(t) = 0, which will be the so ution 
for a system without symmetry breaking. The other solution, relevant for spinodal inflation, is given by the condition 



(v m+mt))). - E ^ (^) n v (2n+1) (j> a 



) = 



(4.1) 



where the subscript a indicates the asymptotic solution. Unlike the sum rule in the large N limit (T™ , this condition 
does not correspond to massless field modes Rather, for a bounded potential of the form ( |2.1[ ) with a definite 
minimum at finite values for the effective mass of particle modes appearing in ( 2.1 2| ) will be positive. Asymptotically, 
we therefore expect the field modes to be redshifted away due to expansion of the universe such that the quantity 
(ijj 2 ) becomes small and may be neglected in (4.1). In this case, the expression for the stationary solution becomes 
simply 



v'(<t> a ) = , 



and the effective mass of the field modes is 



M lff 



v"(<p a ) 



(4.2) 



(4.3) 



These, of course, are just the classical vacuum solutions in t he s ymmet ry b roken phase. 

The task is t o co nnect the early time solutions of Eqs. (3.2) and ( p\4| ) and the late time s olut ions provided by 
Eqs. ( [i2j ) and (4.3). It is convenient to introduce the effective condensate field defined as in (3.7). Doing so, and 
neglecting the exponentially suppressed gradient term in the expression for the energy density (2.23), leads to the 
following equations of motion for the mean field <fi, the condensate field a, and the scale factor a: 



oo 1 

* + 3 ~ & + E -^^v 2 ^ 2 ) (0) = o , 

a ' 2 n n\ 

n— 



(4.4) 
(4.5) 
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1 °° _2n 

2 ^ 2 n n\ 



(4.6) 



Remarkably, these equations are just those one would derive from a classical system of two homogeneous fields with 
the potential 



oc 

n^) = E^ 2n ^ (2n) w- 



(4.7) 



An important property of this potential is that, for interacting fields, it does not depend symmetrically upon <f> and 
a. This means that, in contrast to the large N case fL6| |, it is not possible to combine the mean and condensate fields 
into a single effective classical inflaton. The fields 4> and a have different dynamics and properties. 

The initial condition <r(to) — Hq/2tt results in two distinct regimes. In the classical regime, characterized by 
<P(to) 3* Hq/2tt the dynamics is dominated by the evolution of 4> and is effectively independent of a. This follows 
from the fact that a never grows to be particularly large before <j> reaches its classical minimum. However, in the 
fluctuation dominated regime where 4>{to) -C Hq/2t:, a has a significant effect on the evolution of <fi and dramatically 
modifies the overall inflationary dynamics from naive expectations. 

Before moving on to specific examples, let us examine some general features of the fluctuation dominated regime. 
At intermediate times the dynamics of the reassembled fields will be primarily dictated by the equation of motion 
of (j, due to the fact that a 3> <f> and that both fields contribute to the equations of motion in a similar way. The 
dynamics will therefore approach a quasi-equilibrium regime for which the third term in Eq. (4.5) becomes small: 



(4.8) 
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where the comma represents the partial derivative. Returning to the equation of motion for the mode functions (2.12), 
we discover that this condition corresponds to effectively massless quanta. What we are seeing is the rapid departure 
of the inflaton field from the unstable regime for which its field quanta have a negative mass squared into a physical 
regime of massless quanta. This flattening of the unphysical regime into a form for which particles become well defined 
is reminiscent of the famous Maxwell Construction describing the convexity of the thermodynamical equilibrium free 
energy of a system. The behavior descri bed here is the out of equilibrium analog of such a constructi on fl5|| . 

We note, however, that this condition (4.8) does not correspond to the late time classical solutions (4.1) we expect; 
it is only a quasi-equilibrium for which the mean field 4> continues to evolve toward the true minimum, as is clear 
from Eq. ( gj ). 

To understand the implications of this behavior better, it is useful to work with definite models. We now turn to 
concrete examples which will allow us to follow the dynamics in detail. 



A. New Inflation 



The simplest spinodal model of inflation is new inflation for which the scalar potential ( |2.lD is truncated at quartic 
order with the cosmological constant contribution K = 3/z 4 /2A: 



3/U 4 1 2 2 
— if^ 

2A T 



A_ 

24 



-$ 4 



This model is of particular importance due to its renormalizability. 
The reassembled equations of motion are 

4> + 3-0- /i 2 + V + \o 2 ^ = , 
a 6 2 

CI o A o A , n 

a + 3-cr - ix u + -a 6 + -dfa = , 
a 2 2 



(4.9) 



(4.10) 
(4.11) 
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(4.12) 



The reassembled two field potential for this model is therefore 



3a 4 



60 V) 



(4.13) 



A plot of this two dimensional potential is shown in Fig. |2|. It is characterized by a local maximum at <f> — 0, a = 0, 
a saddle point at = 0, a = fx^/2/X, and minima at (j) = ±/Ua/6/A, <t = 0. 

As we discussed above, there are two dynamical regimes determined by the initial value of <f>. In the classical regime, 
a never plays a significant role and the field (f> acts as an ordinary classical field in a A</> 4 potential ([0]). However, 
in the fluctuation dominated regime, the evolution proceeds first through an inflationary phase with energy density 
contribution given by the tree level potential, but then enters a second regime for which the condition ([T^) is satisfied. 

For this potential, we have in the quasi-equilibrium regime the sum rule 



A 
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, 



(4.14) 



which we recognize as the conditi on f or massless quanta. This expression may then be substituted back into the 
equation of motion for the <p field (4.4), where we find 



A 



-4? = 



(4.15) 



Here we see that the potential energy contribution to the equation of motion for (f> appears at the cubic order. 
The field cf> therefore evolves as a field with an effective squared mass given by — A</> 2 /3. As this is typically much 
smaller (in absolute value) than /x 2 , what we observe is that the potential becomes flattened as a consequence of the 
non-perturbative growth of fluctuations. 







FIG. 2. The two dimensional potential which may be used to describe the full non-linear field dynamics in a A<3? 4 new inflation 
phase transition in terms of effectively classical fields (j> and a. Axes are scaled such that the true minima of the potential occur 
at ((> = ±1. 



We can also use the condition ( 4.14 ) to determine the value of the potential energy ( 4.13 ) during this phase. We 
find simply 
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(4.16) 



as is consistent with the effective mass for 4>. We see that the growth of fluctuations has produced an effective 
dynamics corresponding to a very flat potential for <p with a cosmological constant contribution to the energy density 
/i 4 /A which is 2/3 of the value appearing in the original potential. We therefore arrive at a second stage of inflation 
with an expansion rate related to the original stage by a factor of -\/2/3. 

We note, however, that there is a continued instability in <j) with the consequence that eventually the condition 



> <g = 2 M 2 /A , 



(4.17) 



is met and it becomes impossible to satisfy the sum rule (4.14). At this point, the second inflationary stage ends, 
fluctuation field a decays away and the field 4> evolves to its classical minimum. 



the 



B. Natural Inflation 



The natural inflation potential is derived from the vacuum manifold of a complex scalar field spontaneously broken 
at the Planck scale and with explicit symmetry breaking at the Grand Unified scale. It may be written in the form 

V(<$>) = A 4 [l + cos($//)] , (4.18) 



where A ~ Mqut and / ~ Mpi are constants. Expansion of the cosine reveals that this potential is of the form (2.1) 
with /i 2 = A 2 // and A = A 4 // 4 . 
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The reassembled equations of motion become: 

A 4 

/ 
A 4 



3H 4> — sin ( — ] exp [ — ] = 



(j- 

2/2 



a + 3H& - cos I - j exp \^jj2 ) a = 



(j 



(4.19) 
(4.20) 
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We recognize these equations as those of a two homogeneous classical scalar fields with potential 



(4.21) 



V((f>,a) = A 4 1 + cos £ exp - — 



2/ 2 



(4.22) 



We provide a plot of this two field potential in Fig. pi As expected, for any integer j, we have degenerate maxima 
at (f> — 2jirf, a — and degenerate minima at <j> = (2j + l)irf, a = 0. The feature to notice, however, is that the 
potential quickly becomes very flat as a/f becomes greater than 1. 




Again, we concern ourselves with the fluctuation dominated regime. We find that the sum rule (4.8) results in the 
condition 



A 4 



f 2 
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(4.23) 



which for < irf becomes satisfied as a grows large. The effective mass term for <f> in (4.19) becomes exponentially 
suppressed as well, again indicating that the growth of fluctuations of the field results in a flattening of the potential. 
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The potential itself clearly goes to the value A 4 which is precisely half of the value of the original cosmological constant 
contribution. This state of affairs will continue until <j> > <j> a = nf, at which point the fluctuations represented by a 
become massive and begin to decay away. 

As in the A<& 4 model, we expect two stages of inflation, this time with the expansion rate of the second stage 
reduced from that of the first by the factor \/2. 



V. METRIC PERTURBATIONS 



The possible link between inflationary expansion of the universe at Grand Unified energy scales and the observation 
of fluctuations in the Cosmic Microwave Background (CMB) temperature of order 10 -5 is remarkable |32| . The fact 
that we are now in the process of observing the details of this temperature spectrum through a number of ground-, 
air-, and space-based experiments, presents an amazing opportunity for probing details of what the universe was like 
at times inaccessible through any other means |33|. 

In order to take full advantage of this opportunity, however, it is important that we are careful in connecting 
the observations to the theoretical models of the dynamical processes of the Early Universe which may have led 
to the CMB anisotropics. We have already presented a detailed analysis of the cosmological dynamics of scalar 
fields undergoing inflationary phase transitions. Here, we provide computations of the primordial spectra of density 
perturbations which arise from such phase transitions. 

By computing the spectrum of scalar perturbations, we explicitly show that coherent effects due to large wavelength 
fluctuations of the inflaton field could have significant impact on observational features of the temperature spectrum 
of CMB anisotropics. Particular features of interest include scales on which there is deviation from a flat Harrison- 
Zcl'dovich primordial power spectrum with power increasing as one moves to smaller scales. Such a spectrum, referred 
to as having a blue tilt (corresponding to a scalar tilt parameter n s > 1), was previously thought only to be possible in 
more complicated, multi-field models of inflation. Here, we find that it is indeed possible to produce such a spectrum 
on the length scales relevant to CMB anisotropy observations from the simplest of single field models. 

Another significant feature of the spectra produced in phase transitions is the possible decrease in the amplitude 
of the primordial perturbations on the scales of interest relative to the predictions of an analysis assuming a classical 
evolution for the inflationary field. This may somewhat alleviate the fine tuning problem, although we find in practice 
that this effect is relatively minor, allowing perhaps a dimensionless quartic coupling for the field one or two orders 
of magnitude larger than previously thought. Since this coupling is typically thought to be restricted to be less than 
10~ 14 , we still require the inflationary models to be highly fine tuned. 

A final feature of the spectra is their dependence on the precise initial state of the inflationary field in the region 
corresponding to today's observational universe. It is found that a universe which began the inflationary phase in a 
strongly classical state will show none of the exotic features described here. This strongly contrasts with the case of 
initial states in which quantum fluctuations of the field are of the same order as the classical field value. In this latter 
case, the primordial spectrum may depend quite strongly on precisely how 'quantum' the particular initial state is. 

We begin with a computation of the amplitude of scalar and tensor metric perturbations resulting from spinodal 
inflation. For the specific models discussed above, we provide details of the resulting power spectrum as a function of 
scale. We also include plots of the amplitude and tilt in the power spectrum as a function of the initial state of the 
scalar field for a choice of scale consistent with those observed by the Cosmic Background Explorer |52|. Finally, we 
provide examples of the Ci temperature anisotropy spectra JTs| ] that result from particularly interesting examples of 
spinodal inflationary effects to allow comparison with other theoretical plots of the Ci spectrum and with the observed 
spectrum. 

Note that in what follows, we use the normalizations of the scalar and tensor amplitudes of Ref. 



VI. THE PRIMORDIAL SPECTRUM 



Our st artin g point for the computation of primordial density perturbations is the expression for the average energy 
density ( 2.23 ). To compute the average fluctuation we define 



1/2 



where the variation of the energy density yields the expression 

K) = (^ + ^ + i(wv)H^. 



(6.1) 
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It is convenient to introduce the Fourier variable 



fc 3/2 



defined such that 



<v> 2 



dk, 



In Fourier space, we find the following expression for the average energy density fluctuation on a scale corresponding 
to a wavenumber k: 

5V((f> + 1/)) 



((ST°) (5T$)) k = |V^| 2 (</> 2 + (^ 2 >) + + Ml) ( (tf + i 



dip 



\ipk 



(6.2) 



where terms proportional to 1/a 2 have been neglected. Normally during inflation, the third term on the right hand 
side of (S.2) dominates, which is equivalent to the statement that the usual slow roll conditions are satisfied during 



the evolution. 

Using the above expression, the density contrast at first horizon crossing is 



|*fc= 



1 
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(6.3) 



Here, we note that we have used the important relation between the gauge invariant Bardeen variable representing 
the scalar metric perturbation and the energy density fluctuation at horizon crossing to allow us to write down this 
apparently simple expression. As the density perturbations are adiabatic (recall that there is only one scalar field 
with one set of fluctuations) , the super- horizon evolution of the metric perturbation up to the second horizon crossing 
is specified by the conservation of a single parameter £ as defined in Ref . E3] . Following the procedure of Ref. jOl , 
we arrive at the expression for the density contrast at mode re-entry in terms of the scaled mode functions ipk- 



<i> 2 + <^ 2 >) iv^i 2 + (a + (Mt + mi) + ( (^^) 2 ) m 



25 



(6.4) 



where each quantity is evaluated at the time when the co rres ponding mode first crosses the horizon, i.e. when k = aH. 

In what follows, we will assume that the third term in (6.2) d omi nates over the first two terms, as this is an excellent 
approximation in the models in which we are interested. Eq. (6.4) becomes simply: 
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(6.5) 



This can be simplified further by recognizing that at first horizon crossing the quantity \ipk\ 2 is given approximately 
by (H/2n) 2 fjij . This may be seen directly from the asymptotic solutions of the mode functions for large momenta. 
As the expansion is rapid, these asymptotic solutions remain approximately valid all the way to k ~ aH . Combining 
this with the (semi-classical) Friedmann equation and using the inflationary condition (V) ^> <fi 2 + (tp 2 ), we reach the 
result 
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If we were to make the further assumption that the fluctuations are always small with i/j <C <f> and ip <fi, then we 
would arrive at the standard slow roll expression for the density contrast: 
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We will continue, however, with the more general form ( |6.6| ). 

The computation of the tilt parameter n s — 1 is straightforward, given (6.6). We define 



(6.7) 



d(]n\8 h (k)\) 



(6.8) 



k=aH 



As is common practice, we could rewrite this expression in terms of partial derivatives of the field variables. However, 
since we have dependence both upon 4> and the fluctuations ip, such a proce dure results in a complicated expression 
which is not particularly instructive. We will therefore satisfy ourselves with (6.8), which is used directly to compute 
the tilt parameter in the numerical examples. 



A. A<3> 4 Spinodal Inflation 



Using (BJ3), we now compute the expression for the primordial spectrum of scalar metric perturbations specific to 
the A$ 4 model of spinodal inflation. In terms of the reassembled variables cf> and a, the result is: 
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Again, we point out that all expressions are to be evaluated when the given scale k crosses the horizon. 



B. Natural Inflation 



For natural inflation, the relevant expression is: 
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(6.10) 



C. Gravitational Wave Perturbations 



It is also of interest to examine the spectrum of gravitational wave perturbations resulting from spinodal inflation. 
As such perturbations do not directly interact with the inflaton field, they may be related directly to the expansion 
rate during any inflationary stages. During such regimes, the amplitude of gravitational waves is simply Jl8[ |: 



where Hk is the value of the expansion rate when the scale k first crosses the horizon, fc = aHk- 

As we have discussed, spinodal inflation may involve two distinct inflationary stages. The relevant amplitude of 
the gravitational waves is therefore typically determined by which stage is in effect 60-50 e-folds before the end of 
inflation, when the relevant length scales exit the horizon. However, the transition period between inflationary stages 
may also be relevant. 

The spectrum in the transition period depends upon the details of the transition, but clearly must smoothly 
interpolate between the two major regimes. The important factor is whether the characteristic accelerated expansion 
continues to hold throughout the transition or if the transition includes a short period of deceleration. In the former 
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case, the transition is smooth and straightforward, following Eq.(6.11) throughout, while in the latter case, the 
amplitude of the modes crossing the horizon during the transition has an oscillatory nature |l8f . 

As we will see, the former, smooth transition without any oscillation of the gravitational wave amplitude is typical 
of spinodal inflation, with the result that the spectrum may have at most a single feature indicative of the transition 
from the initial inflationary phase to the spinodal phase. 



D. Notes on metric perturbations 



It is worth taking a moment to examine the significance of the expressions for the amplitudes of the metric per- 
turbations in some detail. The first feature to notice is that the scalar amplitude depends directly upon both the 
average field value <j> and the typical fluctuation represented by a. As we have seen, the end of inflation depends upon 
the evolution of cf>, occurring when </> reaches the spinodal value 4> s . When the influence of the field fluctuations are 
neglected, the value of cf> 60 e-folds before the end of inflation, </>6o, which we will take to be the largest scale measured 
by COBE, is a well defined quantity and corresponds to a well defined scalar amplitude \Sqo \ independent of the initial 
conditions. 

However, we see here that the dynamics of the field fluctuations a may influence the evolution of <f>. The net result 
is that the precise values of 060? ^60: and therefore \5qq\ occuring when the relevant mode crosses the horizon does 
generally depend upon the initial conditions for the inflaton. 

This leads to an extremely important result: The spectrum of primordial metric perturbations resulting from a given 
model of inflation depends upon the initial state of the inflaton field. 

A second significant feature coming as a result of these more general expressions for the perturbation amplitude is 
that |<5fi(fc)| is not necessarily a monotonically increasing function of length scale. In contrast to the case in which the 
influence of the field fluctuations are neglected, there may be periods of time during which |<5/j(fc)| increases as ever 
shorter length scales cross the horizon. 

This leads to a result previously thought to be impossible in such simple, single field inflationary models: Inflation 
may result in a blue spectral tilt in the primordial spectrum of scalar perturbations. 



E. The Ci Spectrum 



As a final note in this section, we recall that none of the parameters of the primordial spectrum written down to this 
point are directly observable. Rather, this primordial spectrum becomes processed by the evolution of the late time 
universe, after corresponding length scales have re-entered the particle horizon. The most important set of observable 
parameters is the C/ spectrum of the Cosmic Microwave Background map of temperature variations. 

Standard techniques of computing the C/'s from the primordial parameters assume a scale independent tilt n s , which 
does not generally apply here. It is, of course, possible to compute the Ci spectrum for a more general primordial 
spectrum, but to do so is far beyond the scope of the present work. We therefore use a standard approximation 
which relates the tilted spectrum cf* 1 ^ to the tilted scalar perturbation amplitude \Sh(k) tdt \ evaluated on the scale 
k = lH^/2 where is the inverse Hubble radius today relative to the respective quantities where a flat spectrum 
(n s = 1) is assumed |3q|: 



^m) _ n { f iat) \s h (k = iH»/2y M 

\8 h (k = lH*/2)f 1 ^ 



Assuming the spectrum C^ lat ' used here is properly COBE normalized, we require that the flat and tilted primordial 
amplitude match for the low I scales corresponding to COBE. 

This will allow us to present an approximate plot of the Ci spectrum resulting from spinodal effects which may be 
compared to standard Ci spectra. 



VII. CONCRETE EXAMPLES AND RESULTS 



We now put all the pieces of the previous sections together in numerical simulations of the full dynamical equations 
of motion of the spinodal system. A note is in order regarding these simulations. In certain cases it was impractical 
to run numerical simulations using the full field theory equations of motion. These include plots showing quantities 
as a function of the initial condition for <fi as well as the plot of the scalar tilt as a function of the initial Hubble 
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parameter. In these cases, we turned to the classical 2-field models described above, using the approximate initial 
condition a(to) = Hq/2tt. All such figures specify in the caption that they result from the classical models. All other 
figures were produced from full field theory simulations. 



A. New Inflation 



We begin with the A<E> 4 system, which we first examined in this context in Ref. fllffl . We show the dynamics of the 
mean field <fi, the fluctuation (-0 2 ) 1 / 2 , each of which is scaled by the factor / = fiy/6/X, and the expansion rate H for 
each of three examples, corresponding to (a) </>(in) 3> H a /2ir, (b) 4>(to) ~ H /2tt, and (c) </>(£n) H /2ir. 

In the first of these, (a), shown in Fig. |^, the evolution proceeds exactly as would be expected from a purely 
classical analysis of the dynamics. The two-point fluctuation (ip 2 ) remains small and does not have a noticeable effect 
on the evolution of 0, which in turn simply follows the contour of the tree-level effective potential V(<fi). As (ip 2 ) 
remains small, the expressions for the amplitude of primordial metric perturbations ([TJ) reduces to the usual slow 
roll expression (3.7) and we arrive at the standard results for |<5h(fc)|, n s — 1, and \5 g \. These quantities are depicted 
in Fig. H as a function of N, the number of e-folds before the end of inflation at which the corresponding length scale 
crosses the horizon. 
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FIG. 4. The mean field (j>(t)/f, the fluctuation < i> 2 (t) > 1/2 //, and the Hubble parameter H(t)/fi vs. t for the A$ 4 model 
with 4>(t Q ) = 5.0i?o/27r, 4>(t ) = 0, Ho = 2/Lt, A/8tt 2 = 1CT 16 , and / = 
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FIG. 5. The scalar amplitude 8h, the scalar tilt n s — 1, and the tensor amplitude 8 g vs. N corresponding to the evolution of 
Fig. | 
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Example (b), Fig. g, is an intermediate example for which the fluctuation becomes large for a short time, inducing 
a spinodal phase during the period of evolution for which the length scales of relevance to CMB observations cross the 
horizon. This example is of further importance because it depicts the case in which the initial classical value of the 
inflaton field is of the same order as its initial vacuum fluctuations. As is clearly shown, the growth of the two-point 
fluctuation has a significant influence on the evolution of the mean field 4>, resulting as well in a modified behavior 
for the expansion. 
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mean field <j>(t)/f, the fluctuation < ip 2 (t) > 1,/2 //, and the Hubble parameter H(t)/fA vs. t for the A$ 4 model 
4H /2tv, 4>(to) = 0,H = 2/Li, A/8tt 2 = 1CT 16 , and / = fijofx. 



FIG. 6. The 
with 4>(to) = 0. 

The result for the metric perturbations are provided in Fig. [?], where we see some remarkable features. First, we 
notice that the amplitude of |5/i(A;)| is reduced and that its shape is significantly changed by the spinodal dynamics. 
The significance of the shape is further emphasized by the scalar tilt n s — 1. Here, we see that for some scales of 
observational relevance between 60 and 50 e-folds before the end of inflation, the scalar tilt parameter n s — 1 is 
positive. This corresponds to a blue tilt in the power spectrum and was, until now, considered to be unattainable in 
inflation models consisting of only a single scalar field. Finally, the tensor amplitude clearly follows the behavior of 
the expansion rate H as expected. 
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FIG. 7. The scalar amplitude Sh, the scalar tilt n s — 1, and the tensor amplitude S g vs. N corresponding to the evolution of 
Fig. | 
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To show why this is the case, we plot in Fig. || the equation of state w = p/p as a function of time. We see that in 
the transition region between the initial and spinodal phases of inflation, there is little departure from the de Sitter 
equation of state p = —p and that, in particular, p < —p/3 through to the end of the second inflationary stage. This 
means that the condition for accelerated expansion d > is satisfied the whole time, resulting in the simple behavior 
of the tensor amplitude. 
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FIG. 8. The equation of state w(t) = p(t)/e(t) vs. t for the parameters corresponding to the evolution of Fig. [| The 
expanded region shows the transition between the initial and spinodal phases of inflation. 

Case (c) , Fig. ^ shows the other extreme case when the initial classical value of the field is significantly smaller than 
the effective quantum fluctuation. Here, (ip 2 ) reaches the spinodal while <p remains small, resulting in a very long 
spinodal phase in which the mean field <f> evolves slowly along the spinodal line. The results for the metric parameters 
are provided in Fig. fLOl 




FIG. 9. The mean field (j){t)/f, the fluctuation < ij) 2 (t) > 1/2 //, and the Hubble parameter H(t)/p vs. t for the A$ 4 model 
with cf>(to) = 0.05#o/27r, 4>{t ) = 0, H = 2p,, A/8tt 2 = 10" 16 , and / = H-J&fX. 
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FIG. 10. The scalar amplitude 5h, the scalar tilt n s 
of Fig. § 
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1, and the tensor amplitude 8 g vs. N corresponding to the evolution 



To add completeness to this picture, it is useful to plot the scalar amplitude and tilt as a function of the log of the 
initial value </>(io) over the range of the differing regimes. This is done for the particular scale crossing the horizon 60 
e- folds before the end of inflation (Fig. fTD) . This figure provides a nice summary of our primary results. 




27tc|>(t )/H 

FIG. 11. The scalar amplitude ($60 and the scalar tilt neo — 1 corresponding to the scale crossing the horizon 60 e-folds before 
the end of inflation vs. 2-K(f)(to) / ft for A$ 4 new inflation using the classical two field model with <j>(to) = 0, A/8-7T 2 = 10~ 16 and 
several values of -Ho- 
rn the classical regime </>(to) 3> Hq/2tt, the results are independent of initial conditions, but as 0(io) becomes of 
order Ho/2n there is a distinct transition regime in which the amplitude drops by up to a couple of orders of magnitude 
(possibly alleviating the fine tuning problem of inflation to a minor degree), while there is a spike for which the tilt 
becomes positive. While this region appears rather restricted in the parameter space of initial values for 0, we remind 
the reader that we plot here on a log scale the ratio of <j>(to) to h/2-k so that the region in question effectively covers 
the entire range for which 4>(to) and Hq/2tt are of the same order of magnitude. 

Finally, we see distinct regions for which there is a long spinodal phase. In this case, (ip 2 ) has a significant effect 
long before the relevant scales cross the horizon. The mean field evolves slowly along the spinodal line to the spinodal 
point. As the sum rule (4.E) is in place during this phase, the evolution of 4> in the latter stages is always the same, 
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and we therefore expect that \8eo\ and uqq — 1 become constant throughout this regime, as is indeed the case. 

It is worth mentioning how these results depend upon the other parameters of the model, such as the magnitude 
of the quartic coupling A and the ratio of the initial expansion rate to the mass scale Hq/ n. 

Increasing the value of A while keeping Hq / fi fixed only acts to modify the value of the vacuum expectation value 
and the spinodal point to be smaller as both of these quantities are proportional to 1/ VX This simply reduces the 
amount of time the fields evolve, but does not change any of the features described above, with the exception that 
the scalar amplitude |<5^(/c)| scales, as usual, as yX. 

Modification of the ratio Hq//i (keeping A fixed), on the other hand has a more complex influence on the behavior. 
The results for the scalar amplitude and tilt as a function of initial conditions for three values of Hq/ are shown in 



Figs. 11 



Here, we see that increasing the parameter has the effect of reducing the amplitudes of the features in the 
transition region. 

Note that in the case of H /[i = 1 the approximations we have used to compute the metric perturbation spectrum 
break down in the transition region, where we see that the scalar tilt spikes to a very large value, and therefore this 
computation is unreliable for that small range of initial conditions. However, as long as the transition period does 
not occur within the last 60 e-folds of inflation, and in particular for the spinodal region for </>(io) ^ Hq/2tt, the slow 
roll conditions are satisfied throughout the relevant period of evolution such that the portions of this plot away from 
the transition region are reliable. 

Finally, we mention a somewhat unexpected feature of the results for the spinodal regime with small initial (j)(to)- 
In Fig. we plot the scalar tilt parameter n s — 1 as a function of the value of the initial Hubble parameter Hq / /i 
and compare it to the result in the classical regime for which n s — 1 = [i 2 /3Hq. The result is remarkable. Over a 
large range of initial expansion rates, there is very little change in the value of the scalar tilt. It seems that the tilt 
due to a long spinodal phase of inflation is relatively independent of the parameters of the model, with an empirical 
prediction over the range of tested parameters of 0.94 < n s < 0.98. While extension to even larger values of Hq/ 
may bring the upper limit toward 1.0, it seems likely that the lower limit will remain solid. The consequence is that 
a future measurement which restricts the tilt to n s < 0.94 would seem to rule out a spinodal phase of inflation. 
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FIG. 12. The scalar tilt neo — 1 corresponding to the scale crossing the horizon 60 e-folds before the end of inflation vs. Ho/fi 
for A$ 4 new inflation using the classical two field model and spinodal initial conditions, 4>(to) "C Hq/2-k. 



B. natural inflation 



Due to its popularity, we examine one other spinodal inflation model in detail: natural inflation. The qualitative 
behavior is in many ways quite similar to that of the A$ 4 model described above, so we will focus on the distinguishing 
feat ures . The primary difference occurs during the spinodal phase. In the A$ 4 model, the spinodal condition of 
Eq. ( 4^ ) simply results in the lowest order contribution to the effective mass squared of the true zero mode <fi being of 
order — A0 2 , thus slowing the evolution of <fi. In natural inflation, however, the spinodal condition for natural inflation 
yields an effective mass which is exponentially suppressed by the growth of the fluctuations. Hence, the evolution 
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of the zero mode can come to practically a standstill. The net result is that the spinodal phase in natural inflation 
models is significantly longer than a corresponding A$ 4 model. 

We provide results for two examples corresponding to the intermediate case with cj)(to) ~ Hn/2ir shown in Figs. |i~3| 
and [l4| and the fully fluctuation dominated case of <f>(to) -C Hq/2-k as depicted in Figs. [l5| and [To] . We see the 
same features that appeared in the previous model, with the obvious difference that the spinodal phase of Fig. [l5| is 
extremely long even for (/>(io) n ot much smaller than Hq/2t:. 




FIG. 13. The mean field <j>{t)j f, the fluctuation < ip 2 (t) > 1//2 //, and the Hubble parameter H(t)/fi vs. t for the natural 
inflation model with 0(t o ) = H /2ir, 4>{t ) = 0, H = 2fi, and A 4 /4?r 2 / 4 = 10" 16 . 
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FIG. 14. The scalar amplitude 5h, the scalar tilt n a — 1, and the tensor amplitude <5 g vs. N corresponding to the evolution 
of Fig. |l| 
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FIG. 15. The mean field <j>(t)/f, the fluctuation < ip (t) > 1 ' 2 //, and the Hubble parameter H(t)/n vs. t for the natural 



inflation model with 0(t o ) = 0.5H /2n, <f>(t ) = 0, H = and A 4 /4tt 2 / 4 
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FIG. 16. The scalar amplitude Sh, the scalar tilt n a — 1, and the tensor amplitude S g vs. N corresponding to the evolution 
of Fig. 



It is worth examining this point further by plotting the number of e-folds of inflation as a function of the initial 
condition on <j). This is shown in FigjlT] where we see a dramatic dependence which contrasts very sharply with the 
logarithmic dependence of the same quantity on initial conditions in the classical regime. 
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FIG. 17. The total number of e-folds of inflation N vs. 2iKf>(to) / Ho for natural inflation using the classical two field model 
with <j>(t ) = 0, A 4 /4vr 2 / 4 = lCT 16 and H = 2A 2 //- 

Once again, it is enlightening to examine the dependence of the quantities |<5eo| and n^o — 1 as a function of the 
initial state. These are depicted in Figs. |l8| where we see again the clear transition regime around 4>(to) ~ Hq/2tt 
separating the classical and fluctuation dominated regimes. 
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FIG. 18. 

field model with <j>(to) = 0, A 4 /4tt 2 / 



1 vs. 2ir<f)(to) / Ho for natural inflation using the classical two 
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and H = 2A //. 



C. Example Ci Spectrum 



Finally, we present an example Ci spectrum resulting from spinodal inflation, shown in Fig. |l9|. We plot the 
tilted Ci spectrum corresponding to the natural inflation case of Fig. [l^ as well as a standard flat Ci spectrum for 
comparison. The 'flat' spectrum was produced using the program CMBFAST version 2.4.1 |36| using the defaults of 
a standard Cold Dark Matter cosmology with a tiltless adiabatic p ower spectrum n s = 1, while the 'curved' spectrum 
was reproduced from the flat one using the approximate relation (6.12). 

The primary feature of interest is that while the spinodal inflation spectrum is shifted downward for the modes 
100 < I < 1000, the spectrum approaches its flat counterpart for high /. This is indicative of the shift of the spectrum 
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from red to blue over the range of observable scales. 
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FIG. 19. The standard CDM C; spectrum of temperature anisotropies vs. 
spectrum corresponding to the evolution of Fig. [L3|. 



I for a flat n s = 1 spectrum and the curved 



VIII. CONCLUSIONS 



Whether one studies QCD, ferromagnetism, or ordinary gases, non-linear, long wavelength effects are seen to have a 
dramatic impact on the properties of the physical system. Without close contact between the observation and theory 
of these systems, there is good reason to doubt whether the respective phenomena of confinement, ferromagnetism, 
or the liquid-gas phase transition would be as well understood as they are today. But that is precisely the challenge 
before us if we wish to understand the dynamics of phase transitions in the very early stages of the universe. 

If we expect to be able to proceed, then it will be necessary to return to those systems that we believe we understand, 
say in condensed matter physics, and examine the techniques, approximations, and concepts which have led to 
consistent and accurate theoretical descriptions. Looking in any statistical physics text we immediately notice that 
there are a few general concepts and techniques which have been found to be particularly rewarding for a variety of 
different systems. 

One of these is mean field theory, where complex details of a system are replaced by simple averages. It is in many 
ways a very naive approach and often provides quantitative results which are only roughly correct. But the real power 
of mean field theory lies not only in its quantitative predictions, rather in the qualitative pictures it paints, allowing 
us, for example, to explain the liquid-gas phase transition in terms of a system as simple as that given by the van der 
Waals equation of state. 

Another powerful concept is the convexity of the thcrmodynamical free energy function for any equilibrium system 
and the simple rules for writing down such a function provided by the Maxwell construction. The story being told 
by these concepts is that theoretical energy curves with concave portions describe 'unphysical' states, and that the 
dynamics of the corresponding system will act to transform any such state into a physical state described by a purely 
convex energy function. 

This report is an attempt to understand the system of an inflationary phase transition in the context of these 
powerful concepts. The initial state is unphysical, described by a concave effective potential, the analog of the 
equilibrium free energy, with quanta corresponding to imaginary mass states. Such a system must decay into physical 
states, which we see from the exponential growth of long wavelength fluctuations. This decay continues until a 
non-perturbative state is reached for which the corresponding quanta are physical with a strictly real or zero mass. 

There are two possible ways to reach such a state. The first is simply to provide a significant bias to the state such 
that the mean value of the field (i.e. the order parameter) reaches its true vacuum value where the field quanta are 
well defined. 

The second way, appropriate for systems with small order parameter, is to allow the system to phase separate into 
domains for which the field has either positive or negative value. Rather than the order parameter moving along the 
potential energy diagram (see Fig. [I]), the field drops down into the center of the diagram to the spinodal line. At 
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the spinodal line, the field quanta are massless and physical, and this state of affairs may be relatively long lasting, 
ending only when one phase becomes so much more prevalent than the other that the system may relax into a definite 
vacuum state throughout the system. 

We therefore have arrived at a consistent physical picture of inflationary phase transitions based upon mean field 
theory Already, at this naive level, we have seen new phenomena which impact not only the evolution of the inflaton 
field but also have important implications for the interpretation of recent and soon to come observational data. 

As a final note, we emphasize that the techniques used here are really only a first - or, rather, a second - approxi- 
mation to a very complicated system of interaction between an unstable scalar field and gravity at very high energies. 
There are a number of possible avenues which might be taken to improve upon these results for the dynamics and, in 
particular, for the predictions of observational quantities. 

One direction is to move beyond mean field theory. Gravitationally, this means doing something more sophisticated 
than semi-classical gravity. There has been recent work in this regard within the context of perturbation theory up 
to two loop order p7[ , and interest in this area has grown somewhat due to the possibility of new phenomena in the 
context of preheating [p8[ . However, the non-perturbative dynamics of the scalar field studied here corresponds to 
non-perturbative departures of the gravitational dynamics from that of a purely classical background field so that 
techniques based upon perturbative expansions are of little help. 

In terms of the scalar field dynamics, one possible avenue that has received attention is the 1/N expansion of the 
O(N) vector model, which includes contributions beyond mean field theory at next to leading order (this approximation 
is also promising as it might be consistently implemented for gravity as well) |39[ . Another alternative approach is to 
use variational methods to compute the dynamics of the system pp| ] , a technique which might also be combined with 
the 1/N expansion. However, significant hurdles remain before either of these techniques will be implementable for 
interesting field theory problems. 

The primary alternative to these semi-analytic approaches is to look to lattice simulations. As the interesting 
phenomena occur out of equilibrium, fully quantum simulations are ruled out and our only alternative is to examine 
real time simulations of a classical scalar coupled to classical gravity. This approach has led, in particular, to 
improvements in our understanding of preheating dynamics However, extending these methods to interesting 

problems in inflationary dynamics will be a challenge requiring a fully general relativistic lattice coupled to the inflaton 
field, a way to properly deal with an exponentially changing scale, and a method of introducing classical fluctuations 
consistent with the classicalization of quantum, sub-horizon field modes as they cross the horizon. 

The challenges are great, but we should be encouraged both by the improved physical picture presented here and 
by the incredible fact that inflationary theory is on the verge of becoming an observational science. 
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